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Abstract 

A recently developed technique for the determination of the density of partition 
function zeroes using data coming from finite-size systems is extended to deal with 
cases where the zeroes are not restricted to a curve in the complex plane and/or come 
in degenerate sets. The efficacy of the approach is demonstrated by application to 
a number of models for which these features are manifest and the zeroes are readily 
calculable. 



1 Introduction 



The study of phase transitions is central to statistical mechanics. Of primary interest is 
the determination of the location, the order and the strength of the transitions. While 
only systems of infinite extent display such phenomena, these are not directly accessible 
to the non-perturbative computational approach, which is restricted to a finite number 
of degrees of freedom. There are, however, well-established techniques for the extraction 
of information from numerical studies of finite systems, and prominent amongst them is 
finite-size scaling (FSS). 

The FSS hypothesis is based on the premise that the only relevant scales are the corre- 
lation length of the infinite-size system and the linear extent of its finite-size counterpart 
£Q . A modification, in which the correlation length of the finite system replaces its actual 
size, extends the validity of the hypothesis to the upper critical dimension Tradi- 
tional techniques to determine phase transition strength from finite-size data involve the 
application of FSS to thermodynamic quantities or to the lowest lying partition function 
zeroes 

However, a full understanding of the properties of the infinite-size system requires 
knowledge of the density of zeroes too. While it has long been expected that extraction 
of this quantity from finite-size systems would be a lucrative source of information, a 
technique to do so proved elusive 0j. The source of the difficulties is that it involves 
reconstruction of a continuous density function from a discrete data set, or sets, as the 
density of zeroes for a finite system is essentially a set of delta functions. Recent consider- 
ations have bypassed these difficulties jSE]. Rather than focusing on the density of zeroes 
itself, one determines the integrated density of zeroes. The robustness and efficiency of 
this approach was demonstrated in E] and the method favourably compared to other 
techniques in j7j. 

In these previous analyses, the distribution of zeroes had two special properties. These 
are (i) the zeroes dominating critical or pseudocritical behaviour lie on a curve called the 
singular line, which impacts onto the real axis at the transition point and (ii) these zeroes 
are simple zeroes (zeroes of order one). While these two properties are common to the 
bulk of models in statistical physics and in lattice field theory, they are by no means 
generic and the question of the generality of the technique presented in jSJ |S] therefore 
arises. 

The purpose of this paper is to extend the method presented in jSJ E] to deal with 
situations where the above two properties do not hold. Instead, the method developed 
here in Sec. 2 assumes the zeroes to be distributed across a two-dimensional region in the 
complex plane and/or to occur in degenerate sets. Such distributions of zeroes have been 
observed in various models of statistical physics and lattice field theory in two dimen- 
sions. The models we address in Sec. 3 are (a) the Ising model on a square lattice (using 
Brascamp-Kunz boundary conditions) with anisotropic couplings, (b) the Ising model on 
a bathroom-tile lattice, and (c) the case of free Wilson fermions in two dimensions. While 
all of these models are in the same two-dimensional Ising universality class, their detailed 
distributions of zeroes are quite different and provide a sufficiently wide sample to test 
the improved density-of-zeroes approach to the detection and characterization of phase 
transitions. Finally, Sec. 4 contains our conclusions. 



2 Zeroes and their Densities 



All of the information on a thermodynamical system in equilibrium is encoded in the 
zeroes of the appropriate partition function. Indeed, for a system of finite size, when the 
partition function, Zl, can be written in as a polynomial in an appropriate function, z, 
of temperature, field or of a coupling parameter, we may write 

Z L (z) = A(z)l[(z - Zj (L)) , (2.1) 

3 

where L denotes the linear extent of the system, j labels the zeroes, and A(z) is a smooth 
non-vanishing function which plays no crucial role in the sequel and is henceforth dis- 
carded. 

In numerical approaches to critical phenomena, FSS of the zeroes, Zj(L) (with j fixed 
- typically to j = 1, which labels the zero nearest the transition point), is used to 
determine properties of phase transitions. A summary of the status of some of these 
calculations is given in [3J. On the other hand, attempts have also been made to gain 
a deeper understanding of some more tractable models analytically [HIE]- Where these 
attempts have involved zeroes of the partition function, it is clear that much information 
is contained in their density. The technique developed in [3] is essentially a convergence 
of these two approaches, and we summarize it here for convenience. 

2.1 Simple Zeroes on a Singular Line 

The reduced free energy is obtained from (|2.1|) as 

f L (z) = ±lnZ L (z) = ^j:in(z-z j (L)) , (2.2) 

3 

having discarded the regular contribution coming from A(z). Here V represents the 
volume of the system. In independent series of publications, Abe |8j and Suzuki [H] 
assumed that the zeroes fall on a singular line in the complex plane, parameterized by 
z = z c + r exp (i(f)) , where z c is the transition point. In this case, a necessary and sufficient 
condition to achieve the correct scaling behaviour for the specific heat is that the density of 
zeroes along the singular line in the thermodynamic limit behave as 5 , o( r ) oc r 1 ^", where 
a is the usual critical exponent of the specific heat. Integrating, gives the cumulative 
density of zeroes in the infinite- volume limit, 

Goo(r)ocr 2 ^ . (2.3) 

In the finite-volume case, the density of zeroes is a string of delta functions, and, 

9l{t) = ^Y,Kr ~ r {L)) , (2.4) 

3 

where the j th zero is given by zj(L) = z c + rj(L) exp(i<j)). Integrating this along the 
singular line leads to the following expression for the cumulative density of zeroes IE] : 



(2.5) 



The two central observations of jSJ |H] were, firstly, that equating the infinite-volume 
density formula ()2.3|) to its finite- volume counterpart ()2.5|) is sufficient (with hyperscaling) 
to recover standard FSS expressions (indeed, FSS, traditionally the consequence of a 
hypothesis, emerges quite naturally from this approach), and, secondly, that (|2.5|) is 
a sensible definition of the cumulative density of zeroes in the finite case. With this 
definition, the strength of transitions may be directly measured by fitting to the ansatz 

G(r) = ai r a2 + a 3 . (2.6) 

In particular, a non-zero value of 03 indicates a definite phase. When 03 vanishes, a 
transition of first order is indicated if 02 ~ 1, while a value of a?, larger than 1 indicates 
a second-order transition with strength 

a = 2 - a 2 . (2.7) 

In [3] and [0] this method was tested by application to a number of models in statis- 
tical physics and in lattice field theory. In all of these models, the locus of zeroes is 
one-dimensional, with a singular line impacting on the real axis at the transition point. 
Furthermore, all zeroes for finite lattices were simple zeroes (with no degeneracies). The 
question now arises as to how the technique translates to more general distributions of 
zeroes. 



2.2 General Distribution of Zeroes 

Departures from such smooth linear sets of zeroes were first observed for models on 
hierarchical and anisotropic two-dimensional lattices, for which there can exist a two- 
dimensional distribution (area) of zeroes ^OlEJE]- Since then, a host of systems have 
been discovered with this feature El [13 EH]- A common characteristic of all such 
two-dimensional distributions of zeroes is that the only physically relevant point at which 
they cross the real axis, in the thermodynamic limit, is that which corresponds to the 
phase transition. It is, however, possible that the zeroes cross the real or imaginary axis 
at unphysical points. These points may be associated with new universality classes. 

Stephenson jTTj has shown that the density of zeroes for such two-dimensional distri- 
butions in the infinite- volume limit is 

9oo(x,y)=y 1 ^- m f(^ ) j , (2.8) 

where (x, y) give the location of zeroes in the complex plane, with the critical point 
as the origin. Here m is a new type of exponent which is related to the shape of the 
two-dimensional distribution [T7j . 

Integrating out the x-dependence in (12 .8|) yields [17] 

f X2 1- 
9oo(y) = g QO (x,y)dx ocy a , (2.9) 

J XI 

where x\ and x 2 mark the extremities of the distribution of zeroes at a distance y from 
the x axis in the complex plane. Integrating again, to determine the cumulative density 
of zeroes at the point r in the y-direction, yields 



Goo(r) oc r 2 



(2.10) 
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Figure 1: Schematic plot of cumulative density of zeroes as defined by (j2.12|) . In this 
example, where the volume, V, is fixed, t\ = r 2 are 2-fold degenerate, while r 3 = . . . = r 6 
are 4- fold degenerate. 

an expression identical to (|2.3|) . The strength of the transition, as measured by a, can 
therefore be determined by similar methods to those previously used. However, rather 
than counting the zeroes along the singular line, one now counts them up to a line y = r 
within the two-dimensional complex domain they inhabit. 

The second new feature we wish to accommodate is the existence of degeneracies in 
the set of zeroes. If a number of zeroes coincide, Gl, as defined in (|2.5jh is multivalued 
and is no longer a proper function. A more appropriate density function is determined as 
follows. Suppose, in general, that Zj = zj +1 = . . . = ^ +n ._ 1 are n^-fold degenerate. By a 
glance at Fig. [T] it is easy to convince oneself that the densities to the left and right of an 
actual zero are given by 



The density at the n^-fold degenerate zero, rj : is again sensibly defined as an average: 



This is the most general formula for extracting the density of any distribution of zeroes 
and deals with two-dimensional spreads and degeneracies. Fitting this quantity to the 
form ()2.6|) yields the strength of a second-order transition through ()2.7|) . As in [S] and 
jH] , the criteria for a good fit are good data collapse in L (or V) and j near the transition 
point and 03 be compatible with zero. 

The error estimates appropriate to this modified density analysis may be determined 
from a procedure adapted from [5] and which we now elucidate. In the present case, where 
zeroes may be degenerate, the monotone nature of the cumulative density function means 
that the actual value of (^(r^L)) cannot deviate from (j2.12|) by more than ±rij/2V (see 
Fig. P). The quantitative difference between this starting point and that in j3] is that this 
deviation is not constant in this case. 

Let G° hs (L) represent the data point coming from the size-L lattice and corresponding 
to the j th zero, which is n^-fold degenerate. Assign an initial error (Jj(L) = cr^rij/V to 




(2.11) 




(2.12) 



this data point, where cr ar b is arbitrary. With these errors, the appropriate goodness-of-fit 
is given by 



Xi = E lGf ' (L l S (L)? = £ J^IGTW - GTW? , (2.13) 

L,j L,j a a.rb n j 

where the expected density value, G^ xp (L), comes from the model ()2.6|) . Minimizing x\ 
yields the parameters a, in (|2.6|) with associated errors denoted daf h . 

Assume, now, the actual error associated with each data point is, in fact, arij/V. The 
corresponding chi-squared may be written [3] 



Xl = ^fxl ■ (2.14) 

If the model fits well, xl/^dof should be close to unity, where iVdof is the number of degrees 
of freedom. The error assigned to each point now becomes 

cr 2 = al bX l/N doi = xl/N doi , (2.15) 

having chosen <j ar b to be unity. Moreover, the actual errors associated with the parameters 
cii are (with cr ar b = 1) 

da, = —da™ h = ada™ h . (2.16) 

Oarb 

Just as in [S], this approach prohibits an independent goodness-of-fit test. 

In summary, the proceedure is to let <Jj(L) = fij/V and minimise x\ i n (|2.13|) to find 
a, and da^ Th . The best estimates for the errors are, then, da, = \Jxx/N& i daf h . 

Note that standard FSS is for fixed-index zeroes and gives that the distance of a zero 
from the critical point is 

r j(L) ~ L~ x l v . (2.17) 

Typically one uses the imaginary part of the zero, Imzj, for the distance Tj in a traditional 
FSS analysis. The real part of the lowest zero may be considered as a pseudocritical point. 
Its scaling is characterised by the so-called shift exponent, A, and 

Re^i(L) - z c ~ L- x , (2.18) 

where z c marks the critical point. Usually A coincides with 1/z/, but this is not always 
the case and the actual value of the shift exponent depends on the lattice topology. For 
a summary of some recent results concerning the finite-size shifting of the pseudocritical 
point in the Ising case in two dimensions, see 



3 Testing the Method on Various (Ising) Models 

We take three two-dimensional models for which the zeroes are easily calculated. In each 
case the real, physical, critical point is in the Ising universality class, with strength of 
transition given by a = (corresponding to a logarithmic divergence in the specific heat). 



3.1 Square Lattice Ising Model with Anisotropic Couplings 



The task of analytically solving the Ising model in two dimensions for finite-size systems is 
greatly ameliorated by the usage of Brascamp-Kunz boundary conditions jTS| EI| > where 
for an M x 2N lattice, the spins in the left boundary row at x = are fixed to + and in the 

right boundary row at x = M + 1 to the alternating sequence H 1 — . . ., whereas in the 

^/-direction periodic boundary conditions are assumed. In the general case of anisotropic 
couplings - Ji along the x- and J2 along the y-direction, with arbitrary ratio R = J2/J1 
- the partition function takes the form [20] 

M N 

Z Mt2N = 2 2NM J] [J [cosh(2/3) cosh(2i2/3) - sinh(2/3) cos(0;) - sinh(2ii!/3) cos(^)] , 
i=i j=i 

(3.1) 

where 4>i = in I (M+l), 0j = (2j — l)ir/2N, and (3 = Ji/fc^T. Recall that for fully periodic 
boundary conditions, the analogue of (|3.1|) consists of a sum of four product terms j^I] 
which is much more cumbersome to analyze for the zeroes. 

For isotropic couplings with R = 1 the term in square brackets of (|3.1|) simplifies to 
1 - 2£ sinh(2/5) + sinh 2 (2/5), with -1 < £ = (cos + cos 6j) /2 < 1. It immediately follows 
that the complex zeroes can be parametrized exactly as sinh(2/5) = £ ± %\J\ — £ 2 , i.e., 
that they are distributed on the unit circle in the complex sinh(2/?)-plane. 

For the anisotropic model with R = 3, each factor in (|3.1|) can be rewritten as a 
fourth-order polynomial in w = 2 sinh(2/3) to give 



J M,2N 



M N 

nn[ 



w 



+ 5w 2 + 4 — 2w cos(0j) — (6w + 2w 3 ) cos( 



(3.2) 



The zeroes of (13. 2 j) are also easily determined numerically, but it is not possible to 
parametrize them by a single variable, implying that they are distributed across a two- 
dimensional region rather than on a one-dimensional curve as in the isotropic case. In 
Fig. Elthis two-dimensional distribution of zeroes is shown for the case R = 3 and a square 
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Figure 2: The partition function zeroes in the complex w = 2sinh(2/?) plane for the 
anisotropic (J2 = 3Ji) L = M = 2N = 40 Ising model with Brascamp-Kunz boundary 
conditions. 
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Figure 3: The distribution of zeroes near w — 1 for the anisotropic Ising model with 
anisotropy ratio R = 3 subject to Brascamp-Kunz boundary conditions for L = 40 — 140 
and j = 1 (x), j = 2 (+), j = 3 (*), j = 4 (o), j = 5 (□), j = 6 (o), j = 7 (.), j = 8 (*). 



lattice of size L = M = 2N = 40. 

The zeroes impact onto the real axis at the point w = 1 and the critical behaviour is 
expected to be dominated by the zeroes close by. The zeroes in this case are all simple 
zeroes (no degeneracies), so it should be noted that this case is essentially a test of the 
applicability of the method to the situation of a two-dimensional distribution of zeroes in 
the complex plane rather than a test of how the method copes with varying degeneracies. 

The cumulative-density distribution for this set of zeroes is plotted in Fig. [31 A three- 
parameter fit to (J2.6j) for the first 8 zeroes for lattices of size L = 40, 60, 80, 100, 120, and 
140 gives as = 0.000002(15), indicating the presence of a transition. With 03 set to zero, 
a two-parameter fit then yields = 2.016(32), close to the expected value of 2 (which 
corresponds to a = 0). 

A closer inspection of Fig. [3] shows that the j = 1 zeroes (denoted by the symbol 
x) are slightly misaligned with respect to the higher-index zeroes. We have therefore 
repeated the fit restricted to j = 2 — 8, which yields G(r) = 0.088(7)r 2 - 008 ^ 33 \ so that 
a = —0.008(33). This is nice confirmation that the technique works when the distribution 
of (non-degenerate) zeroes is two-dimensional. 

Standard FSS applied to fixed-index zeroes using (|2.17|) yields the expected result, 
v = 1 ^H]- Similarly, the shift exponent in (J2.18|) is found to be A = 2. Thus A is not 
coincident with 1/u. This contrasts with the case of the Ising model in two dimensions 
with toroidal boundary conditions [22] but matches results using topologies with a trivial 
fundamental homotopy group [23]. 

To understand these numerical results we return to the finite lattice expansion of (J3.2J) 
and look at the finite-size scaling of the lowest zero wi, which is given by the roots of the 
factor in (|3.2j) with % = j = 1 on an M x 2N lattice. For an infinite lattice the expression 
factorizes to give (4 + w 2 )(l — w) 2 and we see the points where the distribution pinches 
down as the roots at w — 1 (and at w — ±2z). For a finite square lattice (L = M = 2N) 
we can expand around the root at 1 in powers of 1/L to find 

iri 7r(2i + 57r) 

Separating the real and imaginary parts yields lmwi(L) ~ L~ x and Rewi(L) — w c ~ L~ 2 . 




Figure 4: The partition function zeroes in the complex w = 2sinh(2/?) plane for the 
anisotropic (J2 = 2J\) L = M = 2N = 40 Ising model with Brascamp-Kunz boundary 
conditions. 



For comparison, we present a similar analysis with anisotropy ratio R = 2, for which 
the zeroes are plotted in Fig. 0] While the overall shape of the distribution is the same 
as in the R = 3 case of Fig. |2l its detailed structure is different. In this case the density 
analysis reveals a 3 = —0.000 01(2), and a subsequent two-parameter fit to the first 8 
zeroes for lattices of size L = 40-140 yields a 2 = 2.009(30), i.e. a = -0.009(30). The 
corresponding data is displayed in Fig. 
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Figure 5: The distribution of zeroes near w = 1 for the Ising model with anisotropy ratio 
R = 2 subject to Brascamp-Kunz boundary conditions for L = 40 — 140 and j = 1 (x), 
3 = 2 (+), j = 3 (*), j = 4 (o), j = 5 (□), j = 6 (o), j = 7 (.), j = 8 (★). 



3.2 Bathroom- Tile Lattice 

It is also possible to obtain two-dimensional distributions of zeroes for two-dimensional 
Ising models with isotropic couplings, one example being the Ising model on a bathroom- 
tile lattice [inj. This is the (4 • 8 2 ) lattice depicted in Fig. Eland which is dual to the 
Union Jack lattice. 



Figure 6: The bathroom tile lattice. 



The continuum form of the (reduced) free energy on the bathroom-tile lattice is given by 
/ = | + Iln(l + „) (3.4) 

11 2 \n[A{u) + B{u){cos{e 1 ) + cos{e 2 )) + C{u)cos{e 1 )cos{e 2 )] , 



-7T J — 7T (27r) 2 

where u = exp(— 2f3) and 

A(u) = (1 + w 2 ) 2 (l -Au+ 10m 2 -4m 3 + m 4 ) , 

B(u) = 2u(l -n) 3 (l + u)(l + u 2 ) , (3.5) 

C(u) = -Au 2 {l-uf . 



This system is described in detail in [13J. The zeroes of the partition function were 
calculated from the finite lattice discretization of one of the terms in the partition function 
for periodic boundary conditions leading to (|3.4j) . namely 

m n , 

nn U(u)+b(u) 

, 1 a — 1 



=1 s=l 



/2r-l\ (Is -1 

cos — — — + COS 



V M J \ N 



/Or — 1 N /2s — 1 N 1 X / 2 

+ C(m)cos(— — )cos(— — ) \ . (3.6) 



\ M J \ N 

In principle the full partition function is a sum of four 1 such terms, differing in the 
arguments of the cosines which correspond to the four possible choices of (anti)periodic 
boundary conditions for the two species of fermions in the continuum limit of the model. 
In using (J3.6|) . we are assuming that the scaling behaviour of one of these terms is generic. 
An alternative, which we do not pursue here as we are essentially interested in testing 
the scaling of the cumulative density of zeroes rather than formulating the finite lattice 
models themselves, would be to construct Brascamp-Kunz type boundary conditions for 
the bathroom tile lattice. This would also have the effect of projecting out a (different) 
single product term in the expression for Z. 

1 One of which will vanish at criticality for toroidal topology. 
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Figure 7: The partition function zeroes in the u = exp (—2/?) plane for the bathroom-tile 
Ising model (J3IEJ) with L = M = N = 40. Here AFM, PM, FM and Oi indicate the 
anti-ferromagnetic, paramagnetic, ferromagnetic and unphysical phases, respectively. 



The phase diagram for such a system has paramagnetic [PM] , ferromagnetic [FM] and 
anti-ferromagnetic [AFM] phases as well as an unphysical phase which we denote as Oi, 
to adhere to the same notation as ^3]. The zeroes have varying degrees of degeneracy. 
Those for L = M = N = 40 are depicted in Fig. [7|in the complex u = exp (—2/5) plane 
and a blow-up of the region near the ferromagnetic critical point for L = 200 is given in 
Fig. |H| Zeroes in the vicinity of the critical point taper off into a quasi-one-dimensional 
locus, so the bathroom-tile case is a test of the applicability of the method to zeroes of 
varying degeneracies, rather than to a true two-dimensional distribution. 

The physical ferromagnetic critical point is given by u = (1/2) — 2 — \/2 \ = 

0.249 038 4 . . ., corresponding to (3 — 0.695 074 1 . . . In this region, the j = 1 zeroes 
are four-fold degenerate, the j = 5 are eight-fold degenerate, the j = 13 zeroes are 
again four-fold, the j = 17, j = 25 and j = 33 zeroes are each eight-fold degenerate, 
and the j = 41 zeroes are four-fold degenerate. The cumulative density of zeroes near 




Figure 8: The bathroom-tile Ising zeroes near the ferromagnetic critical point u = 
0.249 038 4 . . . for L = M = N = 200. 
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Figure 9: The distribution of zeroes for the bathroom-tile Ising model with L = M = 
N = 40 - 200 and j = 1 - 4 (x), j = 5 - 12 (+), j = 13 - 16 (*), j = 17 - 24 (o), 
j = 25 - 32 (□), j = 33 — 40 (o), j = 41 - 44 (.). 

this ferromagnetic critical point is depicted in Fig. El for L = 40, 70, 100, and 200 with 
j = 1-44 (seven data points for each Li). A three-parameter fit to the form (|2.6|) clearly 
shows that the curve goes through the origin. Indeed, such a fit to the above data 
gives 0,3 = 0.000 000 7(830). Now, setting 0,3 = 0, a two-parameter fit to the data yields 
a 2 = 1.998(18), corresponding to a = 0.002(18), fully consistent with zero, as expected. 

The physical antiferromagnetic critical point is given by u = 4.015 445 4 . . ., near which 
the zeroes again have a one-dimensional locus (as evident in Fig. EJ). The degeneracy 
pattern for the first 44 zeroes is the same as in the above ferromagnetic critical case. A 
three-parameter fit yields 03 = 0.000 01(3), and a two parameter fit to this data gives 
a 2 = 2.03(2). Restricting the fit closer to the origin by using the j = 1 — 16 (3 data points 
for each Li) yields a 2 = 1.9994(163), compatible with a = 0. 

The accumulation point between the ferromagnetic and unphysical regions occurs at 
u = —0.601 231 8 . . . (for which there is no real (3). Here the degeneracy pattern is different 
to those above, with the j '• = 1 zeroes being four-fold degenerate, the j = 5 zeroes eight- 
fold, the j = 13 zeroes again four-fold, the j = 17 and j = 25 zeroes each eight-fold 
degenerate while the j = 33 zeroes are four-fold and the j = 37 zeroes are eight-fold 
degenerate. The density analysis again reveals a transition (a 3 = 0), with a 2 = 2 (e.g., 
the first 24 zeroes for L = 40 - 200 give a 2 = 1.993(12), corresponding to a = 0.007(12)). 

A similar accumulation pattern occurs at the boundary between the antiferromagnetic 
and unphysical Oi phases at u = —1.663 2519, with the corresponding density analysis 
yielding a 2 = 2.0095(123). 

At each of the above four accumulation points, traditional FSS yields v = 1 and A = 2. 

3.3 Wilson Fermions 

The partition function, Zl(k) for a system of free Wilson fermions involves an integral 
over Grassmann variables, which, on completion, leads to the determinant of the Wilson 
matrix, . Here k — 1/ (2m + d) is the hopping parameter, mo is the dimensionless 
bare fermion mass and d is the lattice dimensionality (which is 2 in our case). It is well 
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Figure 10: The partition function zeroes for the L = 50 free Wilson fermions in the 
complex 1/2k plane. 



known that this system exhibits a phase transition at 1/2k = d = 2, where massless 
fermions appear in the continuum limit j21]. This determinant may be expressed as a 
product of eigenvalues, and, for even lattice extent, L, 

Z L (K)=detM^= nilAi 0, (p) , (3.7) 

a=l p 

where 



Ai 0) (p) = ^-Ecosp, + *(-in 



^sin 2 ^ , (3. 



with p^ = 2irp fl /L and where p\ = —(L — l)/2, — (L — 3)/2, . . . , (L — l)/2, while p 2 = 
—L/2, —L/2 + 1, . . . , L/2. These values comply with standard boundary requirements 
for Grassmann variables, namely that they are periodic in the spatial (1-) direction and 
antiperiodic in the temporal (2-) one [24J. 

The complex hopping-parameter zeroes are easily and exactly extracted from the mul- 
tiplicative expression for the partition function (see ^3]) and the zeroes for a system of 
size L = 50 are depicted in Fig. ^] in the complex I /2k plane. 

A special feature of Wilson fermions is the occurence of so-called doubler fermions. 
This means that apart from the physical critical point, which occurs where the zeroes 
accumulate at 1/2k = 2 in the figure, there are lattice artefacts at 1/2k = and at 
I /2k = —2 where further accumulations of zeroes, leading to critical behaviour, occur. 

These Wilson-fermion zeroes clearly form a two-dimensional distribution. They also 
come in degenerate sets, with the first and seventh zeroes being 2-fold degenerate, while 
the third and nineth are 4-fold degenerate. So this system encapsulates both new features 
we seek to address. 

The density plot for the zeroes near the physical transition is given in Fig. ^2 Using 
the first twelve zeroes for lattices of size L = 50, 100, 150, 200, and 250 (four data points 
for each lattice size), a three-parameter fit yields 03 = 0.000 005(29), convincing evidence 
that the density plot indeed goes through the origin. The subsequent two-parameter fit 
yields a 2 = 1.996(11), giving a = 0, as expected. 
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Figure 11: The distribution of zeroes near the physical critical point 1/2/t = 2 for free 
Wilson fermions with L = 50 - 250 and j = 1 - 2 (x), j = 3-6 (+), j = 7-8 (*), 
j = 9 - 12 (o), j = 13 - 16 (□), j = 17 — 18 (o), j = 19 - 23 (.). 

It is worthwhile also applying the method to the artifactual doubler transition at 
1/2k = 0, where the two-dimensional nature of the distribution is more pronounced. 
There, the density data again fall on a universal curve (see Fig. IT2"|) and 03 is determined 
to be 0.000 01(6). A two-parameter fit now yields 02 = 1.996(11), again demonstrating 
that a is zero and the success of the method. Finally, as in the other systems studied 
here, traditional FSS yields v — 1 and A = 2, so in each case the shift exponent does not 
match the inverse of the correlation-length exponent. 
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Figure 12: The distribution of zeroes near the artefactual critical point 1/2k = for free 
Wilson fermions with L = 50 - 250 and j = 1 - 2 (x), j = 3 - 6 (+), j = 7 - 8 (*), 
j = 9-12 (o). 



4 Conclusions 



A recently introduced technique to extract a continuous function, in the form of the 
density of partition function zeroes, from sets of discrete data has been extended to deal 



with the general case where (i) zeroes do not fall on a one-dimensional curve and/or where 
(ii) multiple zeroes may occur. The technique is tested in a variety of models which lie in 
the same universality class as the two-dimensional Ising model and which exhibit various 
combinations of these general features. It is seen to be capable of direct determination 
of the strength of the phase transition, as measured by the critical exponent a. We have 
compared the results obtained from more standard finite-size scaling of the individual 
zeroes and also found good agreement. 

It also perhaps worth highlighting that in this exercise we have found that formulating 
an Ising model with anisotropic couplings and Brascamp-Kunz boundary conditions is 
straightforward and still leads to a simple product form for the finite lattice partition 
function, a very useful property for investigating scaling. Though we have only touched on 
the topic briefly in this paper, the exotic critical points which appear at complex couplings 
in many models are also amenable to our analysis, and we discuss this elsewhere. 
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